Feshbach resonances in ultracold atom-molecule collisions 
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We investigate the presence of Feshbach resonances in ultracold alkali-dialkali reactive collisions. 
Quantum scattering calculations are performed on a new Na3 quartet potential energy surface. An 
analysis of scattering features is performed through a systematic variation of the nonadditive three- 
body interaction potential. Our results should provide useful information for interpreting future 
atom-molecule collision experiments. 
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I. INTRODUCTION 

In recent years, there has been a growing interest in ul- 
tracold molecules [3], particularly in the production and 
properties of the molecules formed from ultracold atomic 
gases HL Photoassociation 0| and Feshbach resonance 
tuning [4j are two main experimental techniques for a co- 
herent production of ultracold molecules from ultracold 
alkali-metal atoms. In 2003 long-lived molecular Bose- 
Einstein condensates were created from weakly bound 
homonuclear lithium and potassium dimers by exploit- 
ing magnetically tunable Feshbach resonances between 
fermionic isotopes [1, @, 0|. 

While Feshbach resonances are always located at the 
highest vibrational manifold of the dimer, photoassoci- 
ation could in principle allow access to low vibrational 
dimer states. In 2005 RbCs molecules were created in 
their ground vibronic state Very recently, several 
different photoassociation schemes for molecular forma- 
tion in the ground vibronic state have also been devel- 
oped [lESEIEiEl. 

For the correct interpretation of the forthcoming ex- 
periments with cold molecular samples it is essential to 
understand the atom-molecule and molecule-molecule in- 
teractions at sub-K temperatures Theoretical re- 
sults have already been published for the homonuclear 
X+X2 ultra-low-energy collisions with X = Li, Na, K (TEI . 
EE EE El, HE Hl|. Isotopically heteronuclear Li+Li 2 
ultra-low-energy collisions have also been studied theo- 
retically [2]], [12, [13, [11, [25[. Collision cross sections have 
been measured experimentally for the Cs+Cs2 ultracold 
inelastic processes [1EII3]- Very recently reactive and in- 
elastic rate constants were measured for Li+Lij at room 
temperature j28|. 

Cold collisions are known to be very sensitive to poten- 
tial energy surfaces [H, Hl[ , and therefore experimental 
information is needed to improve the corresponding the- 
oretical models. In particular, knowledge of low-energy 
resonance patterns often allows different properties of the 
interaction potential to be determined with high accu- 
racy. Such resonances have been studied in great detail 



both theoretically and experimentally in ultracold atomic 
gases; see e.g. (29I [30I l3lj . 

On the other hand, very little is known about atom- 
molecule resonances in ultracold collisions. Similarly to 
atomic scattering the hyperfine-induced resonances could 
in principle exist at very low collision energies. However, 
model calculations have shown that for a general polar- 
ization they will be quenched by inelastic spin-exchange 
transitions forming singlet molecules [32]. 

Alkali-metal dimers on the lowest electronic triplet 
manifold are only stable if they are in a doubly spin- 
polarized state under collisions with doubly spin polar- 
ized atoms (assuming that the relativistic spin interac- 
tions are neglected). Unfortunately, for this specific po- 
larization hyperfine-induced resonances are prevented by 
symmetry (resonances induced by relativistic spin inter- 
actions are still possible). In spite of this, long-lived 
three-atom complexes can in principle exist and give rise 
to resonances. Such reactive resonances have been iden- 
tified in reactive collisions at room temperature [33T ]. 

In this work, we focus on collisions of Na2 molecules 
in the triplet electronic state with ground-state Na 
atoms. For this study a new potential energy surface 
of Na3(l 4 j4 2 ) nas been constructed. The occurrence of 
long-lived three-atom resonances in such collision com- 
plex is demonstrated in the ultracold regime for the first 
time. We also study the dependence of collision cross sec- 
tions on the potential energy surface and we show that at 
least knowledge of two terms in the cross section partial 
wave expansion is needed in order to characterize further 
the three-body potential. 



II. POTENTIAL ENERGY SURFACE 

Ab initio calculations were performed using a single- 
reference restricted open-shell variant [H] of the coupled- 
cluster method [35] with single, double and non-iterative 
triple excitations [RCCSD(T)]. A basis set consisting of 
[12s,12p,5d,2f,lg] basis functions [IB] was used for the 
dimer calculations, and the same basis set without the 



g functions was used for the trimer calculations. Elec- 
trons from the Is orbital on each sodium atom were not 
correlated in the coupled-cluster calculations. The three- 
atom interaction potential was decomposed into a sum of 
pair-wise additive and non-additive contributions 

■imcr 'dimer 
i<j 

(1) 

It has been shown by several authors that in the 
case of alkali-metal trimers the non-additive term 
^23j ^13) is rather large and cannot be neglected 
[33, [33]. Interaction energies were calculated with re- 
spect to the separated-atoms dissociation limit, and the 
full counterpoise correction of Boys and Bernardi (ioT ] 
was employed to compensate for the basis set superpo- 
sition error in both the dimer and trimer calculations. 
All the ab initio calculations were performed using the 
MOLPRO quantum Chemistry package [411 ] . 

47 dimer interaction energies Vdimcr on the a 3 £+ 
manifold were calculated on an irregular grid cover- 
ing the range of interatomic distances from 2.0 A to 
14.0 A. These points were interpolated using the ID 
reciprocal-power reproducing kernel Hilbert space (RP- 
RKHS) method [42]. The interpolation was done with 
respect to r 2 using RP-RKHS parameters m — 2 and 
n = 3). The resulting curve had a minimum at approxi- 
mately r e = 5.194 A, V dlmcr (r e ) = -172.946 cm" 1 , which 
is slightly higher than the previously reported ab initio 
minima r e = 5.192 A, Vdimcr^e) = -177.7 cm -1 [13], 
r e = 5.20 A, FdimcrOe) = -176.17 cm" 1 |37j], and 
r e = 5.214 A, V dimcl (r e ) = -174.025 cm" 1 

Ivanov et al. [44( analyzed experimental data on triplet 
Na2 and derived the accurate position r e = 5.16607 A, 
Vdimerf^e) = —173.64960 cm -1 of the a 3 S+ minimum. 
Therefore our ab initio interaction energies were shifted 
and scaled (shifted by -0.02754 A and scaled by 1.00407) 
so that the minimum of the modified potential energy 
curve coincided with the minimum determined from ex- 
periment. The RP-RKHS interpolation was then re- 
peated using the modified RP-RKHS method (45|. Be- 
yond the last ab initio point, the potential energy was 
then extrapolated to the form 

t r / \ ^6 Cg ClO , , 

Fdhner (r) = - ^ - ^ - . (2) 

The long-range coefficients and Cg were kept fixed 
to the values of 1.561 x 10 3 E h ag and 1.16 x 10 5 E h a%, 
respectively The value of the "free" long-range coef- 
ficient Cio was then determined from the corresponding 
RP-RKHS coefficients [13] to be 1.19 x 10 7 E h a\° ', which 
compares very well with 1.158 x 10 7 Eh aj° from Ref. [i3 ]. 
The resulting potential energy curve supports 16 vibra- 
tional bound states and gives a scattering length of 67.1 
ao, which compares reasonably well (within 10%) with 
published values 65.3 [48], 63.9 f49], and 62.51 \M\- 

356 trimer interaction energies Vt r i mcr on the l 4 A' 2 
manifold were calculated on a regular 3D grid covering 
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FIG. 1: Cuts through the Na3 quartet surface in valence coor- 
dinates. Upper panel: cut for a bond angle of 60°; the global 
minimum of —881 cm -1 is at ri2 = rn — r^i = 4.34 A. Lower 
panel: cut at collinear geometries; the collinear minimum of 
—382 cm -1 is at r\i = rn = 5.06 A. Contours are labeled in 



the range of interatomic distances from 2.5 A to 10.0 
A (geometry configurations were unique up to a permu- 
tation of atoms and satisfied the triangular inequality 
|?'i2 — ri3| < r23 < r\2 + ^13; at linear geometries, where 
^23 =?"i2 + r i3: the distance r23 was permitted to extend 
beyond 10.0 A). 

The grid consisted of 220 C2V points (including 16 D 3 h 
points) and 136 Coo V points (including 16 Dooh points). 
From the trimer interaction energies the counterpoise 
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corrected non-additive energies V3 were extracted using 
Eq. ©. 

The non-additive energy function V3 was represented 
in the same manner as in the case of the spin-polarized 
potassium trimer In order to accommodate 

the geometric dependencies of the long-range multipole 
terms, third-order dipole-dipole-dipole (H3 | and dipole- 
dipole-quadrupole [53j terms were subtracted from the 
non-additive energy V3. Their corresponding long-range 
coefficients C 9 and C\\ were fixed to 1.892 x 10 5 E h [46| 
and 1.46812 x 10 5 E h aj 1 respectively 0]. The lead- 
ing term of the remaining multipole asymptotic expan- 
sion was the fourth-order dipole-dipole-dipole term [551 ] . 
and after a multiplication by a suitable function it was 
prepared for an "isotropic" extrapolation [ToL l5ll |. The 
resulting points were then interpolated, using the fully 
symmetrized 3D RP-RKHS interpolation method [371 ] . 
in each interatomic distance with respect to the reduced 
coordinate p — (-^) and with RP-RKHS parameters 
S = 10.0 A, to = 0, n = 2. 

The three-atom interaction potential Vtrimer for the 
1 4 A' 2 state of Na3 was then reconstructed using Eq. (TJ. 
Its D3/J global minimum -880.9 cm' 1 is at rv2 = ^13 = 
T23 — 4.34 A and Doo/,, saddle point -381.7 cm -1 is at 
ri2 = f'13 = 5.06 A. The minimum of our trimer po- 
tential is approximately 5% deeper than the minima re- 
ported by Higgins et al. [37] and Soldan et al. [38]. Two 
cuts through the surface are shown as contour plots in 
Fig. [TJ for values of the valence angle 60° and 180°. 



III. QUANTUM DYNAMICS 

The scattering observables are obtained by solving the 
time- independent Schrodinger equation for three atoms. 
Quantum dynamical calculation are performed using hy- 
perspherical democratic coordinates. This system of co- 
ordinates comprises of three internal coordinates (two hy- 
perangles and one hyperradius) describing the shape and 
the size of the molecular triangle and three Euler an- 
gles describing the orientation of the molecular plane in 
space. The total wavefunction is expanded on a set of 
hyperspherical basis functions varying with the hyperra- 
dius. The resulting closed coupled equations are solved 
using a log-derivative propa gato r approach. Details on 
the method can be found in [56]. 

The hyperspherical democratic coordinates are espe- 
cially well adapted in describing alkali species reac- 
tions that mainly proceed through an insertion mech- 
anism [la ]. However, the region of large interparticle 
distances, where the system separates into atom and 
molecule, is not efficiently described in hyperspherical co- 
ordinates. Therefore the scattering wavefunction in the 
outer region is computed using Jacobi coordinates. State- 
to-state probability amplitudes are finally extracted by 
matching to the short-range wavefunction obtained with 
the hyperspherical approach (BH ]. 

Atom-atom collisions in the ultracold regime are de- 
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FIG. 2: The Na+Na 2 (t> = 0, j = 0) s-wave scattering length 
as a function of the three-body control parameter A (see text). 
Horizontal lines identify model potentials corresponding to 
a s = — lOOao and 60ao. 



termined by a small number of parameters. In fact, the 
s-wave scattering lengths and the dominant term in the 
long-range multipole potential expansion are sufficient 
in general to predict all near-threshold scattering and 
bound state properties. This feature has allowed many 
systems of experimental interest to be accurately mod- 
eled based on a limited amount of experimental informa- 
tion [H, EO, [HI . Observed energy-dependent cross sec- 
tions and Feshbach or shape resonances often provided 
a key piece of information for determining theoretically 
the scattering lengths and the long-range dispersion co- 
efficients. 

The situation appears to be more complex for atom- 
molecule collisions due to the additional ro-vibrational 
degrees of freedom and anisotropic interactions. As col- 
lisional data may be soon available we begin to study 
here what experimental information might be best suited 
to constrain the collision models. For most alkali sys- 
tems, and for Na2 in particular, the two body poten- 
tial is known with high accuracy from a combination of 
conventional room temperature and ultracold atom spec- 
troscopy [H3|. Therefore, one may expect the three-body 
interaction V3 to represent the largest source of uncer- 
tainty. We assume that its shape is essentially correct 
and following the approach of [If}, we solve the scatter- 
ing problem with a scaled potential A V3. At variance 
with Ref. [16] which considered inelastic scattering we 
discuss here scattering resonances in elastic collisions. 

We focus on molecules in the lowest triplet rovibra- 
tional state, which is collisionally stable under two-body 
collisions with atoms if both colliding partners have max- 
imal spin projection on the quantization axis. We show 
in Fig. [2] the atom-diatom s-wave scattering length a s as 
a function of the three-body control parameter A. Each 
time a three-body bound state crosses the dissociation 
threshold the a s presents a typical divergence, termed a 
zero-energy resonance. One may note that a 0.1-1% po- 
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tential variation (1-10 cm -1 on the potential depth) is 
sufficient for a complete — oo to +00 variation of a 8 . 

In order to investigate the relation between the zero- 
energy quantity a s and finite-energy scattering, we select 
values of the control parameter A corresponding to the 
same value of a s and compare the corresponding energy- 
dependent elastic cross sections. 

We first consider total angular momentum J = 0. For 
j = rotational states this implies an angular momen- 
tum I = 0. Fig. [3] shows the result of the comparison 
for a typical positive value of a s . The partial J = 
cross sections a for a s > show a qualitatively similar 
behavior essentially determined by the value of a s and 
by the long-range Ce coefficient. One can remark the 
well known zero energy limit a — ► Ana 2 . The minimum 
of the cross section corresponds to the scattering phase 
shift going through a multiple of 7T, and in the absence of 
contributions from higher order partial waves would cor- 
respond to a Ramsauer-Townsend minimum in the total 
cross section 1171. 
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FIG. 3: The partial J = Na+Na 2 (^ = 0, j = 0) elastic cross 
section as a function of collision energy for different values of 
the three-body control parameter determined by the same 
value of the atom-diatom scattering length a s = 60ao. The 
inset refers to the typical energy regime of current ultracold 
molecule experiments. 

The cross sections calculated for a s < are more inter- 
esting. A large negative scattering length is associated 
with a bound state turned into a virtual state with energy 
Eo in the continuum. This situation can be conveniently 
described by decomposing the elastic phase shift into a 
background plus a resonance contribution [l7l . 



5 = S hg + 6, 



<5rcs = — arctan 



7/2 
E — Ef 



(3) 



where E r = Eo + r\ is the resonance position, with 7 
and r\ the resonance width and shift, respectively. For 
the first few angular momenta £ low energy scattering 
is determined by the asymptotic behavior 7 ~ E i+1 / 2 : 
S hg ~ E l+1 ' 2 and rj ~ const (H3|. 



There will be a low-energy resonance only if the width 
is sufficiently narrow 7 (E) < E. If the more strict con- 
dition 7 (E) <C E is fulfilled 7 (E) can be replaced with 
the constant quantity 7 r = 7 (E r ) and the decomposition 
in Eq. |(3]) implies that S undergoes a rapid ir variation 
across resonance. Note that because of the 7 ~ E 1 / 2 
threshold law, I = scattering at sufficiently low energy 
always violates the 7 < E condition and no resonant be- 
havior will arise. However, this does not rule out the 
presence of resonances at higher yet very low collision 
energy (see below). 

Resonances can also be analyzed in terms of the 
Wigner time delay [58], 



i.e. the average delay of a scattering event compared 
to free transit in the absence of the potential. In the 
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FIG. 4: The partial J = Na+Na 2 (^ = 0, j = 0) elastic cross 
section as a function of collision energy for different values of 
the three-body control parameter determined by the same 
value of the atom-diatom scattering length a a = — 100ao. 
Cross sections presenting resonant behavior are identified by 
dashed lines. The dotted line denotes the unitarity limit (see 
text). Sample cross sections obtained using different A and 
presenting very similar energy dependence are emphasized in 
the inset. The inset also shows (dash-dotted lines) two virtu- 
ally identical cross sections for a s > extracted from Fig. [3] 



threshold regime using Eq. |(3]) one obtains: 



Q = 



(e - E r y - 

E — E r 

(E - E r ) 2 - 



11 

4 



1 - 



drj 
dE 



2 dSbg 
v dk 



1 d-y 

livdk 

4 



(5) 



where v is the velocity in the relative motion and k 
the relative wave vector. The first term is the usual 
Lorentzian profile arising from exponential decay with 
a ^ correction. Near isolated resonances at high energy 
this is usuallythe dominant contribution to the time de- 
lay; see e.g. [59(|. The resonance shift rj is usually a slowly 
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varying function of energy but in the presence of addi- 
tional scattering features such as shape resonances in the 
background continuum [60]. 

The second term in Eq. |(5]) is the classical time for 
the relative particle to span a distance 2-%r- For 1 = 
elastic scattering <5b g ~ —ka,h g , and this term reduces to 
—2a hg /v corresponding to an attractive (repulsive) char- 
acter of the background potential for negative (positive) 
background scattering lengths a^ g . The third term van- 
ishes for E = E r and gives a correction of dispersive 
shape across resonance. 

Note that in Fig. [4] one can essentially identify two 
classes of curves. The first class (full lines) presents a 
monotonically decreasing behavior towards the first min- 
imum, which corresponds to nonresonant scattering. The 
second class (dashed lines) present peaks at which the 
scattering partial cross section reaches the unitarity limit 
c = |f (dotted line in Fig. [3]) and can in principle be as- 
sociated with a resonant behavior. 




FIG. 5: The Wigner time delay Q (upper panel) and the sin 2 8 
quantity for J = Na+Na2(« = 0, j — 0) elastic collisions and 
selected A values. Only the peak near 400/iK can be classified 
as a resonance (see text). 

We focus on the lowest-energy peak of Fig. [4] and make 
the three-body potential slightly more attractive in or- 
der to further shift this feature towards threshold. Fig. [5] 
shows the sin 2 S and Q quantities for three selected poten- 
tials. The feature near 400^iK can be essentially classified 
as a resonance with j r ~ 0.5E r , As the potential be- 
comes more binding 7 becomes larger than E r (j r = E r 
for E r — 200^K) and the resonant behavior tends to dis- 
appear. As the peak is made to shift closer to threshold 
the time delay is fully dominated by the background con- 
tribution; see upper panel. Also note (lower panel) that 
in all cases the unitarity limit sin 2 S = 1 is attained. 

An additional interesting feature that can be observed 
by inspection of Figs. [30 is the near coincidence of J = 
elastic cross sections computed with different three- 
body potentials in the whole energy range E < 10 mK. 
Sample cross sections illustrating this circumstance for 



both a s > and a s < are put forth in the inset of Fig.[4l 
This shows that knowledge of energy dependent cross 
sections in the regime where only the J = partial wave 
is important is in general not sufficient to determine the 
strength of the three-body nonadditive potential, even if 
its shape were precisely fixed. 

However, model potentials giving equivalent J = 
cross sections will not in general be equivalent if J > 
scattering is explored. For instance Fig. [6] shows the 
J = 1 elastic cross sections for (v = 0,j = 0) molecules. 
Calculations use the same model potentials as the in- 
set in Fig. [4] and can be identified according to the line 
style. Order of magnitude differences are observed for 
,7=1 cross sections in cases where J = cross sec- 
tions are identical. One can conclude that the initial 
characterization of a theoretical model based on purely 
elastic collisions should take into account a sufficiently 
broad energy range for the contribution of J > par- 
tial waves to become observable. In alternative, in cases 
where the method for the production of cold molecules 
allows the initial ro- vibrational state to be controlled, at 
least one experimental inelastic cross section (for some 
(v = 0, j > 0) initial state, for instance) should comple- 
ment the elastic collision data. 
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FIG. 6: The partial J = 1 Na+Na 2 (^ = 0, j = 0) elastic cross 
section calculated with the same set of model potentials used 
in the inset of Fig. [4] Same line style is used for corresponding 
cross sections obtained with the same potential. A resonance 
can be observed in the upper dashed curve. The inset shows 
the resonance width j r as a function of its position E r . 

Please note that one peak is also observed in Fig. [6l 
Analysis of numerical results based on Eq. ((3]) shows that 
7 r ~ 0.5EV, i.e. this feature represents a resonance. The 
inset shows the resonance width as the peak center is 
shifted by making the A parameter vary. Its position 
below the maximum ~ 400/iK of the t = 1 centrifugal 
barrier (we find 7 r = E r for E r ~ 500/iK) suggests that it 
is a shape resonance. However, Feshbach coupling similar 
to the one found for J = collisions is not conclusively 
ruled out. 
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In conclusion, we have presented a new potential en- 
ergy surface for Na3(l 4 A2). We have demonstrated that 
long lived triatomic complexes exist and give rise to res- 
onance effects in reactive collisions even at very low col- 
lision energies. General features to be expected in atom- 
molecule scattering in the ultracold regime have also in- 
vestigated by performing a systematic variation of the 
three-body part of the interaction potential. Knowledge 
of energy-dependent J — elastic cross sections may not 
be sufficient to determine the strength of the nonaddi- 
tive three-body interaction. To this aim, at least one 
additional J > elastic or inelastic cross section needs 
to be experimentally determined. In this work we have 
studied the sensitivity of scattering observables by in- 
troducing a global scaling parameter of the three-body 
interaction. In perspective, as cold collision empirical 
data will become available, it is likely that more complex 



parametrizations of the potential energy surface will need 
to be introduced in order to compare quantitatively the- 
ory and experiments. 
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